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ABSTRACT 

We use the UKIDSS Ultra-Deep Survey to trace the evolution of galaxy clustering 
to z = 3. Using photometric redshifts derived from data covering the wavelength range 
0.3 — 4.5/mi we examine this clustering as a function of absolute K-band luminosity, colour 
and star-formation rate. Comparing the deprojected clustering amplitudes, we find that red 
galaxies are more strongly clustered than blue galaxies out to at least z = 1.5, irrespective of 
rest-frame K-band luminosity. We then construct passive and star-forming samples based on 
stellar age, colour and star-formation histories calculated from the best fitting templates. The 
clustering strength of star-forming galaxies declines steadily from tq ~ 7/i _1 Mpc at z ~ 2 
to tq ~ 3/i~ 1 Mpc at z ~ 0, while passive galaxies have clustering strengths up to a factor 
of two higher. Within the passive and star-forming subsamples, however, we find very little 
dependence of galaxy clustering on K-band luminosity. Galaxy 'passivity' appears to be the 
strongest indicator of clustering strength. We compare these clustering measurements with 
those predicted for dark matter halos and conclude that passive galaxies typically reside in 
halos of mass M Js 1O 13 M0 while luminous star-forming galaxies occupy halos an order of 
magnitude less massive over the range 0.5 < z < 1.5. The decline in the clustering strength of 
star-forming galaxies with decreasing redshift indicates a decline in the hosting halo mass for 
galaxies of a given luminosity. We find evidence for convergence of clustering in star-forming 
and passive galaxies around z ~ 2, which is consistent with this being the epoch at which the 
red sequence of galaxies becomes distinct. 

Key words: Infrared: galaxies - Cosmology: large-scale structure - Galaxies: High Redshift 
- Galaxies: Evolution - Galaxies: Formation. 



1 INTRODUCTION suggest that there are crucial differences in the evolutionary his- 

tories of red and blue galaxies at low redshift. Present day mas- 



The origin of the red sequence of galaxies, and more gen- 



sive, passiv e galaxies are claimed to be largely in place by z=l or 
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(Visvanathan & Sandage 1977), have become increasingly impor- in 7- — J i lUnmk ■ i • i j ■ j c ^ 

: : ^ . I " -| i r i Conselice et al. 2007), implying an extremely rapid period ol star- 

tant questions in extragalactic astronomy. Since Dressier ( 1980) , .„,,.. . . . . , 

1 < ^ •* ■ rnrmQtiiw I \-\t* nmitiAcini nnrl mace tiinprirMic r»t mupr.mQcc Qtirl 



formation. The luminosity and mass functions of lower-mass and 

star-forming galaxies, however, evolve significantly over the same 
found in the denser regions of clusters, evidence has grown to ^ ^^ (Bundy gt ^ 2QQ6 . q j^ 



first identified the trend for elliptical galaxies to be preferentially 



The early formation of st ars in the most ma ssive galaxies has 
E-mail: ppxwhl@nottingham.ac.uk become known as downsizing JCowie et al. 1996) and is in contrast 
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to the well-established hierarchical growth of the underlying dark 
matter distribution. These massive, quiescent galaxies dominate 
present-day groups and clusters and are thought to have formed 
at the highest peaks of the dark matter density field. The expected 
cluste ring of the dark matter is well understood (e.g. IMo & White! 
,2002;), and at any given epoch the real-space clustering ampli- 
tude over mid to large scales is a monotonically increasing func- 
tion of dark matter halo mass. Measuring the clustering properties 
for galaxy samples therefore provides a direct link from the visi- 
ble galaxies to the invisible dark matter in which they are found. 
Furthermore, the growth of dark matter halos over time is also 
well modelled, and so the evolution of the most massive systems 
throughout cosmic history can be traced through the use of cluster- 
ing statistics. 

Over the last few years the clustering technique has been 
applied to complete samples of optically-selected galaxies, sep- 
arated by luminosity and colour. Zehavi et al. (2005) computed 
the angular auto-correlation functions of red and blue volume- 
limited sub- samples of galax ies from the Sloan Digital Sky Sur- 
vey (SPSS; lYork et alj|2000h for separations < 10/i _1 Mpc and 
redshifts, z < 1. They found that the population of galaxies red 
in g — r colour are more strongly clustered than those with bluer 
colours (ro = 5.7 and 3.6/i _1 Mpc for red and blue populations re- 
spectively). Moreover, they show that this difference remains when 
the samples are further broken down into absolute R-band magni- 
tude su b-samples. 

IMeneux et al] J2006t) sp lit the VIMOS-VLT Deep Survey 
(VVDS; lLeFevre et alj|2005l) galaxies i nto red and blue spectral 
types by the method of lZucca et al.l J2006l) and computed their clus- 
tering out to z ~ 1.2. They found that the deprojected clustering 
scale length of the red spectral type galaxies (ro ~ 5ft" 1 Mpc) 
was greater than that of the blue type galaxies (r p ~ 3/i _1 Mpc) 
at all redshifts investigated. ICarlberg et alj J1997|) find sim i lar re - 
sults for their sample at z ~ 0.6. More recently ICoil et al.l (2008) 
used the DEEP2 galaxy redshift survey to address the same ques- 
tion at z = 1. They also found that red-sequence galaxies are more 
strongly clustered (ro ~ 5 and 3h~ Mpc for red and blue sam- 
ples), but additionally showed that those galaxies with intermedi- 
ate colours had an intermediate clustering strength. They computed 
the relative bias between their sub-samples, finding a smoothly in- 
creasing bias with rest-frame (U-B) colour up to the red sequence. 

Recently, two further deep surveys have been used to study 
the clustering of galaxies in the range 0.2 < z < 1.2: t he Canada- 
Franc e-Hawaii Telescope Leg acy Survey (CFHTLS , lllbert et al 



l2006t> and the zCOSMOS field JLillv et alj2007h . McCrack en et al 



(2008) used the CFHTLS, splitting the galaxy sample by template- 
fit spectral type and B-band luminosity. Based on photometric red- 
shifts, they find that the difference in clustering strength between 
galaxies fit by early type and late type templates is roughly constant 
across all redshifts studied. The luminosity dependence at a given 
redshift, however, is minimal. Intriguingly, the early type popula- 
tion at 2 ~ 0.5 shows hints of an inverse dependence of lumi- 
nosity on clustering strength. The authors attribute this behaviour 
to fainter red galaxies being preferentially located in cluster envi- 
ronments. IMeneux et alj ( 120090 meanwhile, use the spectroscopic 
zCOSMOS data, finding a weak dependence on luminosity which 
remains when stellar mass is used instead. 

Each of these studies has been based upon an optical colour 
selection, which limits these studies to z ^ 1. Beyond this red- 
shift optical studies become increasingly biased against passive and 
dusty galaxies and so are only able to probe the clustering of unob- 
scured star-forming galaxies. The relative clustering strengths for 



complete samples of passive and star-forming galaxies are there- 
fore poorly understood beyond z ~ 1. Infrared survey data are not 
biased in this way and in addition, at z < 2, the K-band luminos- 
ity is a reasonable proxy for the stellar mass of a galaxy. Infrared- 
selected samples are therefore vital in the study of high redshift 
galaxy clustering. 

Utilising infrared data, the high redshift progenitors of 
nearby high-mass, passive systems have been suggested as be- 
ing passively-selected BzK galaxies JDaddi et al. 2004: . Kong et alj 
2006; Hartley et al. 2008), ultra-compact z=2 red gal axies (e.g. 
Izirm et all 120071: iToft et alj 120071: ICimatti et al.l 120081) , sub-mm 
galaxies JLillv et alj|l999l : ISwinbank et al. 2006), 24/xm sources 



20081), distant 



selected at z > 1 (e.g. Magliocchetti et al. 

red galaxies (DRG s) jFranx et alj 120031 ; IConselice et al] |2007|; 

Foucaud et jdj 20070. e xtremely red objects (EROs) iRoche et alj 



120021 : iMcCarthv et ah 2004) and Ultra-luminous infrared galax- 
ies (ULIRGS) (e.g. Farrah et al. 2006). Each of these populations 
have been found to cluster very strongly, with rp values in the 



range 7 < rp < 14fc~ 1 Mpc teoche et alj|2 002: Blai n et al 



Farra h et alj|20q6l ; iFoucaud et alj|2007l ; iMagliocchet ti et al 
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In Hartley et al. (2008 ) (hencefor th H08) we used the colour 
selection criteria of lDaddi et alj (|2p04) to investigate the clustering 
properties of passive and star-forming galaxies at 1.4 < z < 2.5. 
We showed that, when limiting the samples to Kab < 23, the pas- 
sive galaxies were significantly more clustered than star-forming 
galaxies. This result was found in angular clustering and also when 
deprojected to find the real-space clustering strength; with cluster- 
ing scale lengths of ~ 15/i _1 Mpc for the passive sample (pBzK) 
and ~ 7/i _1 Mpc for the star-forming sample (sBzK). 

Althou gh all of these d ifferent galaxy populations overlap to 
some extent JLane et a l. 2007), their biases differ and they are sen- 
sitive to different epochs. It is therefore difficult to compare these 
samples with low redshift samples or to study global galaxy evo- 
lution. The analysis in H08 in particular left some open questions. 
We established that passive galaxies are more clustered even at such 
high redshifts, but the passive population is exclusively bright and 
the star-forming galaxies exhibited a strong limiting-magnitude de- 
pendence on clustering strength. It is therefore not clear whether it 
is luminosity or passivity that is more significantly correlated with 
clustering str ength. Even at low-r edshift debate remains over such 
a distinction JNorberg et alj '2002), complicated by the dominance 
of pass ive systems at the b right end of the luminosity function. 

I Williams et alj (2009) used a two-colour rest-frame selection 
to select a quiescent cloud and star-forming track over the range 
1 < z < 2. Using this novel selection, they compute the cluster- 
ing scale lengths of their two samples, finding very similar lengths 
to those of the BzK selections in H08. They also report evidence 
for a mild evolution in clustering with luminosity for the star- 
forming ga laxies by splitting their samples into two luminosity 
bins. ICoil et alj 1.2006) also find, at z = 1, only a modest depen- 
dence for clustering strength to increase with B-band luminosity, 
covering 0.7L* < L < 1.4L*. 

A remaining open question is to determine when the passive 
population is built up. Though w e found in H08 that this build-up 
occurs over 1.4 < z < 2.5, and lWilliams et alj 1 120091 ) find similar 
evidence for their sample at 1 < z < 2, these ranges in redshift are 
very broad. Can we identify a narrow epoch at which the cluster- 
ing strengths of passive and star-forming galaxies of a given lumi- 
nosity are equal? At this epoch passive and star-forming galaxies 
would occupy the same mass of dark matter halos and differences 
in their internal processes could be clearer. In this paper we attempt 
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to disentangle the effects on clustering strength due to passivity 
and luminosity and to further refine the redshift range over which 
the passive population of galaxies is first established. 

Where relevant, we adopt a concordance cosmology in 
our analysis; Q M = 0.3, fi A = 0.7, h= H /100 = 
0.7 kms _1 Mpc _1 and erg = 0.9. All magnitudes are given in the 
Vega system unless otherwise stated. The remainder of the paper 
is organised in the following way: in the following section we pro- 
vide details of the data sets used, derived quantities and method of 
selecting passive galaxies. In section 3 we compute the clustering 
strengths of these samples and then discuss their implications in 
section 4. We conclude in section 5. 



2 DATA SETS, DERIVED QUANTITIES AND SAMPLE 
DEFINITIONS 

The analysis we present in this work is based primarily on the 
third data release (DR3) of the ongoing UKIRT (United King- 
dom Infra-Red Telesco pe) Infrare d Deep Sky Survey, Ultra-Deep 
Survey (UKIDSS UDS. lLawrence et alJ2007h . UKIDSS comprises 
five complementary sub-surveys with the UDS being the deepest, 
targeting a depth of K — 23 over a single, 4-pointin g mosaic of 
the Wide-Field Camera (WFCAM. Icasali et al'|[2007h : an area of 
0.88 x 0.88 degrees. UKIDSS began in the spring of 2005 and the 
UDS has seen several data releases: an early release of the first few 
weeks of data (ED R, iDve et alj|200q). a one-year relea se (DR1, 
IWarren et alj2007h . which we used in lHartlev et alj J2008I) . and the 
more recent three-year data release (DR3). 

The principle improvement from the DR1 to the DR3 is the 
addition of H-band data. Additional data in the K-band also pro- 
vides an incremental improvement in depth. The 5a depths within 
2" diameter apertures for the DR3 in J, H and K-bands are 22.8 
(AB 23.7), 22.1 (AB 23.5) and 21.8 (AB 23.7) respectively. These 
depths make it the deepest single-field near infrared survey of its 
area to date. For details of the stacking procedure, mosaicing, cat- 
alogue extraction and dept h estimation we ref er the reader to Al- 
maini et al. (in prep.) and lFoucaud et al.l 120071) . 

In addition to near infrared data, the field is covered by ex- 
tremely deep optical data in the B, V, R, i' and z'-bands from 
the Subaru-XMM D eep Survey (SXDS, ISekiguchi et"afl 120051 : 
iFurusawa et al.l 12008), Spitzer data reaching 5a depths of 24.2 
and 24.0 (AB) at 3.6/im and 4.5^m respectively from a re- 
cent Spitzer Legacy Program (SpUDS, PLDunlop) and U-band 
data taken with CFHT Megacam. The SXDS utilised Suprime- 
cam on the Subaru telescope, achieving depths of Bab = 
28.4, V AB = 27.8, R AB = 27.7, i' AB = 27.7 and z' AB = 26.7 
( IFurusawa et alj2008l . 3a, 2"). Border regions, areas around bright 
stars and obvious cross-talk artifacts were masked and any sources 
within such regions were discarded. The co-incident area with the 
UDS taking masking into account is 0.63 deg 2 . 

The world-public DR3 catalogue extracted from the K-band 
image is used as the basis for our selection, upon which we impose 
a cut at K — 21.1 (Kab = 23) within 2" diameter apertures. 
This cut was made to minimise photometric errors and the number 
of spurious sources and to ensure a high level of completeness and 
robust photometric redshifts (see below). This limit is fainter than 
My, the characteristic lum inosity of the K-band luminosity func- 
tion dCirasuolo et alj2010l and fig.G}. We are therefore probing the 
properties of 'normal' galaxies out to Z — 3. Corresponding 2" 
magnitudes for each object were extracted directly from the optical, 



J and H-band images at the position of the source after a detailed 
astrometric re-alignment to the UDS K-band image. 

From our combined catalogue we remove stars in the follow- 
ing way: bright (K < 18.1) stars were removed by excluding 
those with K-band aperture magnitudes within 3" and 2" satisfying 
K a — K a > —0.14, where the stars are clearly separated from 
the galaxies. The remaining, fainter star s form a stellar loc us in the 
(B-z')~ 0' - K) plane (as noted bv lDaddi et alj2004l) and can 
easily be removed by t he criterion (z' — K) < 0.3(B — z') — 0.5 
dCirasuolo et al.l 120101) . Saturated stars and the surrounding con- 
taminated regions were carefully masked out during the analysis. 
These cuts left 33923 galaxies from which we made our selections, 
as detailed in the following section. 

The magnitudes obtained from the images were then used to 
determine photometric redshifts (photo-zs) and stellar ages by a \ 2 
minimisation over a large suite of model templates, constructed us- 
ing Bruzual & Chariot (2003) with a Salpeter initial mass function, 
and in cluding a treatment of dust content (see lCirasuolo et al.l2007l 
uOld for a full account). These templates cover nine values of the 
exponential star-formation decay rate (r = 0.1, 0.3, 0.5, 1, 2, 3, 5 
and 10 Gyrs, and an instantaneous burst), and allow stellar ages up 
to the age of the Universe at any given redshift. The dust reddening 
is al lowed to vary betwe en < Ay < 2, following the Calzetti 
law JCalzetti et al.ll2000D. and line of sight neutral Hydrogen ab- 
sorption, calculated from lMadaul (1995J), is also taken into account. 

We remove from our catalogue any source with an unaccept- 
able fit (x 2 > 15, 4.0% of the galaxy sample) as they will likely be 
mis-assigned when we break up our sample by redshift, and the de- 
rived stellar ages are unlikely to be reliable for such objects. Upon 
investigation, the majority of these discarded sources were found 
to be cross-talk artifacts (26%), QSOs (36%) or the minor mem- 
bers of pairs or mergers (23%). The majority of the remaining dis- 
carded objects were of very low surface brightness and many of 
them would not have made it into the volume-limited sub-samples. 
The fraction of otherwise useful objects rejected by this \ 2 > 15 
criterion is therefore < 0.6%. Rest-frame magnitudes and colours 
were computed by integrating the appropriate filter over the best- 
fitting template of each object. In this work we make use of the ab- 
solute rest-frame K-band magnitude and rest-frame (U—B) colour. 



2.1 Sample selection 

In this section we outline our definition of the passive galaxy 
sample. At low redshift, passive galaxies form a well-defined 
red sequence in a diagram of rest-frame U — B colour ver- 
sus absolute magnitude. The red sequence is not exclusively 
composed of passive systems, however, and typically contains 
a 30% contamination from dust y star-forming galaxies (e.g. 
IWolf. Gray. & Meisenheimerl 120051) . At high redshifts (z > 1) 
contaminatio n by dusty star-for ming objects becomes increasingly 
problematic JDaddi et alj|2004l) . In addition, a simple U — B cut 
does not make full use of all the available photometry. For these 
reasons, we define our passive sample with two criteria: (a) the 
galaxy must be 'red' in rest-frame U — B colour and (b) the best- 
fitting galaxy template must imply a very low level of residual star 
formation. These criteria are further detailed below. 

Our rest-frame U — B colours are computed from the best- 
fitting template used to find the photometric redshift. The templates 
cover two possible star- formation histories: either an instantaneous 
burst parameterised by an age, or an exponentially decaying star- 
formation rate parameterised by an age and the e-folding time, r. 
Hence, 
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Figure 1. The distribution of UDS galaxies used in this work in the pho- 
tometric redshift - Mk plane. The dot-dashed lines show the sub-sample 
regions for which we compute clustering lengths, provided they contain 
> 150 galaxies. Also plotted, as the dashed-line curve, is the fit for the 
evolution of the characteristic luminosity, M^. , of the K-band luminosity 
function, as given in Ciras uolo et alj 420100 . 




Figure 2. Colour-luminosity density plot of the full galaxy sample (light, 
red contours) and of those used to define the red sequence (heavy contours). 
The red sequence was defined by fitting a line of the form U — B = a X 
Mk + b to the galaxies which were best fit by a burst template with age 
> 1 Gyr. This fit is shown by the dot-dashed line, while the dotted line is 
the 3<r lower envelope. Galaxies redder than the dotted line are defined as 
'red' for the purposes of our sample selections. 



the burst galaxies described above, the associated best fit (dashed- 
dotted line) and red/blue selection criterion (dotted line). 

The red sequence is know n to evolve with time, with redder 
colours at lower redshifts (e.g. iBrammer et al.1120091) . Accurately 
accounting for this evolution with photometric redshifts is prob- 
lematic. The typical errors on our TJ-B colours are 0.15 magni- 
tudes, which is of the same order as the expected colour evolution. 
For simplicity, we choose not to treat the red sequence in such de- 
tail and instead use a fixed colour selection boundary. We note that 
our 'red' sample includes part of the green valley, and so we are 
unlikely to miss many intrinsically red galaxies at high redshift. 

In order to remove the galaxies that are red due to dust rather 
than age we make use of the star- formation histories obtained from 
the templates. To be selected as a passive galaxy we require the 
galaxy to be red, as defined above, older than lGyr and to have a 
best-fit SFR (from the photo-z template fits) below some cut-off 
value (as given by equationQ] where r = f or a burst model). By 
varying this SFR cut-off criterion we can study how the inclusion 
of galaxies with greater recent star formation rates influences the 
clustering strengths of the sub-samples. 

We require, for our primary passive sample, a highly con- 
servative limit to residual star-formation rate: SFR ^ 0.1% of 
S-F-Rinitiai (in addition to being red in rest-frame U — B and with 
a stellar age > lGyr). We choose this strict limit to minimise the 
contamination by dusty star-forming objects whilst retaining suffi- 
cient objects with which to compute meaningful clustering lengths. 
There are well-known degeneracies between age, metallicity and 
dust content. We cannot completely overcome them, even with the 
multiple bands of deep photometry at our disposal. However, these 
parameters imply a strong break in the galaxy SED and hence only 
the most extreme dusty contaminants will be included. The num- 
ber of objects satisfying these selection criteria is 3991 (33% of the 
'red' sample). The vast majority of passive galaxies selected in this 
way also lie within the 'pERO' colour selection boundaries (see 
section[4]l. 

We choose, as a baseline against which to compare our pas- 
sive sample, an actively star-forming sample with a single criterion: 
SFR > 10% of SFRinMmi- This star-forming sample consists of 
22856 objects of which 91% are 'blue' according to our (U-B) cri- 
terion and 9% are 'red'. We then construct three further passive 
samples by relaxing the SFR criterion used in our primary passive 
sample to include galaxies with SFR < 0.5%, SFR a ^ 1% 
and SFR Sj 10% of S-F-Rinitiai- In each of these three additional 
samples the requirements for red U — B colour and stellar age are 
kept the same as the primary sample. The number of objects in 
each of these expanded passive samples are 6111, 7677 and 10032 
respectively. We will occasionally refer to these samples in the text 
by the SFR limit imposed during selection. 



SFR obs = SF-Rinitiai x e 



-age/-r 



(1) 



where SFR b 3 is the SFR at the time of observation. 

We wish to define our red sequence initially using the most 
conservative and clean sample of passive galaxies drawn from our 
catalogue. We construct this sample to be those galaxies which are 
simultaneously old (age > lGyr) and are fit with a burst template, 
as these values imply a strong 4000A break. We then perform a \ 2 
minimisation to fit an equation of the form U — Bab = ax Mk + 
b to these galaxies. The 'red' sample is then taken to include all 
galaxies within 3<r scatter of this fit (or redder). In Figure|2]we plot 
contours showing where galaxies lie in the (U — B)ab vs. absolute 
K-band magnitude plane. Overlaid in black are the contours for 



3 CLUSTERING PROPERTIES 

The real-space correlation length of a sample of galaxies is inti- 
mately related to the mass of dark matter halos in which they are 
found. The value of this clustering scale-length is therefore of great 
interest in the study of galaxy formation and evolution. In this work 
we determine the angular correlation function of the galaxies as 
they are seen in projection, and then use the redshift distribution 
to deproject to find the real-space scale length (ro). This process is 
detailed below. 

The 2-point angular correlation function, w(9) is defined by 
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the joint probability of finding two galaxies in soli d angle elements 



the joint probability ot hnding two galaxies in solid 
6Q,i and 50,2 at a given separation (Peebles 1980), 

5P = n 2 50 1 60 2 (l + tu(0i 2 )). 



(2) 



To estimate the co rrelation function we use the estimator of 
lLandv&Szalavl(ll993h . 



w(9) = 



DD - 2DR + RR 
RR 



(3) 



where DD, DR and RR are the counts of data-data, data-random 
and random-random pairs respectively at angular separation 8, nor- 
malised by the total number possible. Although this estimator is rel- 
atively robust against systematic errors, there remains a small bias 
due to the finite field size. This bias is corre cted by an integral con- 
straint, constant C. We follow the method of lRoche & Ealesl 119991) 
by using the random-random counts to estimate the size of this 
bias: 



C = 



ZN RR (8)8- C 
XNrr(8) 



(4) 



where the sums extend to the largest separations within the field. 

The angular correlation functions are fit by y 2 minimisation 
with a power law of fixed slope, 8 — —0.8. The fit is made over mid 
to large scales to avoid biases arising from multiple halo occupa- 
tion. The small-scale limit corresponds to 0.3/i~ 1 MptQ where pre- 
vious w(8) measurements of galaxies over 1 < z < 2 begin to di- 
verge from the large scale power law (H08). The large scale cut off 
is at 0.2 degrees where survey area limits our measurements. This 
method is the same as that presented in H08. Correlation function 
uncertainties are determined using bootstrap resampling, which are 
typically 2-3 times the Poisson errors. 

The real space clustering and project ed clus t ering are 
linked by the relativistic Limber equation JLimbeii 1 19541) . If 
the redshift distribution of a sample is known, the Lim- 
ber equation can be inverted and the corr elation l ength , 
rp, can be calculated in a robust manner (Peebles 1980; 
IMagliocchetti & Maddox1ll999l;lRo"che, Dunlop. & Almainii. 2003). 



Following IMagliocchetti & Maddoxl d 19991) and adjusting to take 
into account our known photometric redshift distributions we de- 
termine ro as 



r () 



cA 

HnHry 



(J °° n{z)dzf 



oH-r W °° n{z) 2 x(zy-rP(O ,z)F(z)dz 



l/ 7 



(5) 



of the sub-sample 
meanings as in 



where n(z) is the redshift distribution 
and all other symbol s hav e the sarr 
IMagliocchetti & Maddoxl dl999h . 

Below we present the resulting measurements for ro as a func- 
tion of redshift and K-band luminosity. This is first determined 
for all galaxies and then further subdivided by colour and star- 
formation history. 



3.1 The global evolution of clustering 

Using the method outlined above we computed the galaxy cluster- 
ing scale lengths, ro, for each of the volume limited sub-samples 
shown in Figure[T] provided that each subsample contained at least 
150 galaxies. These measurements are shown in Figure [3] The first 



1 This value was chosen to maximise the fitting range, though we have also 
tried 0.5/i~ 1 Mpc and 1.0/i _1 Mpc, finding consistent results. 
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Figure 3. Upper. The global evolution of clustering in the UDS field. Galax- 
ies are split into the volume-limited sub-samples shown in Figure [T] and 
the real-space correlation lengths are computed (see text). The evolution 
of clustering strength with redshift is mild for most luminosity ranges. 
Lower. The sub-samples are further separated into red and blue in rest- 
frame (U — B) colour and the clustering computed once again. Globally, 
red galaxies are slightly more strongly clustered than blue galaxies. The in- 
crease in ro found for bright galaxies at low-redshift is clearly seen to be 
dominated by red galaxies. 



thing of note in these measurements is the lack of strong evo- 
lution of ro with redshift at a given Mk- Furthermore, there is 
very little luminosity segregation, consistent with the findings of 
ICoil et alJ ([2008), although the most luminous galaxies are slightly 
more strongly clustered at z ~ 1. 



3.2 Colour and star-formation rate dependant clustering 

In much of the literature concerning the evolution of galaxy cluster- 
ing, populations are separated by colour, rather than star-formation 
properties. Rest-frame colours are simpler to compute and there- 
fore allow more straightforward comparisons with previous work. 
In Figure[3]we again compute ro values for the subsamples in Fig- 
ure [T] but this time we also separate red and blue galaxies, using 
our rest-frame (U — B) colours defined in Section l2~T"l (Figure[2"ll. 
At Z < 1, red galaxies are sign ificantly more clu stered than 
blue galaxies, similar to the result of ICoil et al.l ( 120080 at z ~ 1. 
This difference is particularly pronounced for the brightest sub- 
samples, where the strong clustering observed in the full sample is 
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Figure 4. Evolution of the clustering strengths with redshift for our four passive samples with the single comparison star-forming 
sample. With the most conservative passive definition (requiring red rest-frame U — B colour, stellar age > 1 Gyr and SFR ^ 
0.1% S-F.Ri n i t i a i), the passive galaxy sub-samples are significantly more strongly clustered than the star-forming sub-samples for 
z < 1.5. Above 2 = 2 the number of passive galaxies becomes too small to make robust clustering measurements. As the passive 
galaxy selection criterion of SFR (given in each panel) is relaxed, the clustering strengths of star-forming and passive sub-samples 
become more similar. Each point is plotted at the mean redshift of the sub-sample it represents. 



clearly a result of the clustering of red galaxies. Above z — 1.5 the 
clustering strengths of the blue and red subsamples appear to con- 
verge, indicating that this may be the principle epoch at which the 
red sequence is populated. The red population is inhomogeneous, 
however, containing both passive and dusty star-forming objects. 
We therefore use our template fits to separate the passive and star- 
forming galaxies (Section r2.lt . 

We plot the deprojected clustering lengths of our passive and 
star-forming samples in Figure [4] with the SFR limits described in 
section |2~T| shown in each panel. For the most conservative passive 
sample (SFR ^ 0.1%) the clustering strengths of the passively- 
selected galaxies are greater than those of the equivalent star- 
forming galaxies for almost all sub-samples up to z = 1.5. Our 
result, shown in the upper left panel of Figure[4] shows that the clus- 
tering length is more strongly dependent on passivity than on lumi- 
nosity, with very little luminosity segregation in the star-forming 
sample. This difference continues to at least the median redshift of 
BzK-selected galaxies l IHartlev et al. 2008). Above z ~ 1.5, how- 
ever, the respective clustering strengths of passive and star-forming 
galaxies appear to converge. Further survey data will be required to 
confirm this finding, which is discussed in section [4] As we relax 
the SFR requirement to allow galaxies with higher residual star- 
formation rates into the passive selection, we find that the strength 
of clustering of the passive sub-samples decreases, which confirms 
that truly passive galaxies are more strongly clustered among the 
'red' population. 



For illustration, in Figure \5\ we show the angular clustering 
measurements for passive and star-forming galaxies over all sub- 
samples in the two redshift ranges, 1 < z ^ 1.5 and 1.5 < z ^ 2. 
The passive samples use our conservative definition (SFR ^ 
0.1% SF7?initiai) described in section |2~T1 The solid lines are the 
best power-law fits and, in agreement with our earlier measure- 
ments for BzK-selected galaxies (H08), the passive galaxies are 
found to cluster more strongly. 



3.3 Comparison with modelled halo masses 

We use the formalism given in IMo & White! 120021) to find the 
clustering strengths of dark matter halos of given masses as a 
function of redshift. Thes e models gat her to g ether previous work 
jPress & SchechteJ 1 1974 : Bond et al. Il99ll : IMo & White! Il99d 
|jing|l998MSheth. Mo. & To rmen 200 lj) and use the ellipsoidal col- 
lapse model which has been calibrated against dark-matter only N- 
body simulations. Their aim is to provide the abundances and clus- 
tering of dark matter halos as functions of redshift and mass across 
a wide dynamic range. 

In Figure [6] we plot these models for halos of mass 1O 1O M0 
to 10 14 M© together with the results from our conservative pas- 
sive sample (see section |2~H and star-forming sample. At z < 2, 
star-forming galaxies generally occupy halos of mass 1O 12 M0 or 
smaller, while passive galaxies occupy halos of mass 1O 13 M0 to 
5 x 1O 13 M0. The general trend for star- forming galaxies is for 
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Figure 5. Angular clustering measurements for conservative passive galax- 
ies (SFRo Sj 0.1% S-Filinitjai) and star-forming galaxies within the 
redshift ranges 1 < z ^ 1.5 (upper) and 1.5 < z ^ 2 (lower) 
with the star-forming points off set for clarity. Subsamples in luminos- 
ity have been combined for illustrative purposes. Solid lines show the 
best power-law fits, using a constant slope of —0.8, over large-scales: 
0.3 h~ 1 Mpc (Dashedline) — 0.2 deg. In both cases the passive samples 
are more strongly clustered, consistent with the measurements of pBzK and 
sBzK-selected galaxies (e.g. H08). Errors are determined from a bootstrap 
analysis. 



Figure 6. The clustering of our conservative passive sample and star- 
forming sample (symbols have the same meanings as in Figure|4j a re com- 
pared with the clustering predictions of dark matter halos from Mo & White! 
12002). Lines of constant halo mass (in Mq) are shown, with additional 
dashed-dotted lines for 5 X 10 13 A/0 (upper) and 5 X 10 12 Mq (lower). 
Downsizing in the star-forming population is evident with the same lumi- 
nosity galaxies found in less massive halos towards z=0. The same can- 
not be said of the passive galaxies, which typically occupy halos of mass 
M ^ 10 13 M Q . 



galaxies of a given luminosity to occupy halos of greater mass 
at higher redshift. This is a manifestation of halo downsizing 
l lCowie et alJI 19961 : Foucaud et al. 2010, in press). This downsiz- 
ing in host halo mass for star-forming galaxies appears to extend 
beyond z=2. At these very early times bright star-forming galax- 
ies are hosted by massive dark matter halos (M ~ 5 x 10 Mq). 
Presumably these are among the progenitors of low redshift group- 
dominant galaxies. 

The passive galaxies occupy halos up to an order of magni- 
tude more massive than the star-forming galaxies. In contrast to the 
star-forming galaxies, the brightest passive galaxies appear to show 
constant clustering strength over the range 0.5 < z < 1.5. If con- 
firmed by future survey data, this result suggests that bright passive 
galaxies are hosted by similar mass halos across this redshift range. 



4 DISCUSSION 

We have shown that passive galaxies, defined by fitting templates 
over 11 bands from the UV to Spitzers IRAC bands, are more 
strongly clustered than star-forming galaxies at z < 1.5. Clustering 
strength is intimately linked with the minimum mass of dark matter 
halo that can host a galaxy of a given type. Passive galaxies there- 
fore on average occupy more massive dark matter halos than their 
star-forming counterparts. These results indicate that at z < 1.5 
the passivity of a galaxy sample is a strong indicator of the typical 
mass of dark matter halo that host them. 

The K-band luminosity is a reasonable proxy for the stellar 
mass of a galaxy, though galaxies may dim by approximately one 
magnitude after t hey stop forming stars (e.g Lillv & Longair 1984; 
ICowie et alj 19961) . Even taking this dimming into account, the con- 
servative passive sub-samples (SFR ^ 0.1% SFRi n i t i a i) are typ- 
ically more strongly clustered than the star-forming samples for a 
given stellar mass. 



Furthermore, there is an absence of a strong luminosity depen- 
dence on the clustering strength of the star-forming samples. We 
have therefore found that of the two effects that we set out to dis- 
entangle, passivity and luminosity, passivity is the more significant 
indicator of clustering strength at z < 1.5. However, we note that 
we are unable to probe the most extreme high luminosity galaxies 
(Ma- < —26) due to our limited survey area. 

Above z — 1.5 we are unable to significantly distinguish 
star-forming and passive galaxy clustering. The correlation lengths 
of passive and star-forming galaxies apparently converge at this 
epoch. If this behaviour is confirmed it would suggest that the 
epoch z ~ 2 is the epoch in which the passive and star-forming 
samples are first becoming distinct. Hence, it is likely to be the ma- 
jor epoch at which the red sequence is being populated. This finding 
makes the further study of the 1.5 < z < 2.5 range critical and we 
intend to return to and improve upon this work as the UDS data 
push deeper and spectroscopic samples become available. 

In the hierarchical formation of structure the first halos to col- 
lapse are those which eventually merge to form the most massive 
halos at low redshifts. Our result then points towards a time se- 
quence: passive galaxies formed earlier in those first halos while 
those of similar stellar mass, but still forming stars, developed in 
lower mass halos that collapsed later. Galaxy and halo evolution 
is accelerated at the earliest epochs with respect to the present 
day. The galaxies formed at these strongly clustered, high density 
peaks of the matter distribution are therefore likely to become fully 
evolved, passive galaxies more rapidly than the general population. 
Discounting mergers, the stellar mass of a galaxy is limited by the 
available gas reservoir so a passive galaxy could have formed sig- 
nificantly earlier than a star-forming galaxy, but end up with very 
similar stellar masses at the epoch of observation. Their respective 
halos, however, will have built-up mass since they first collapsed 
and so those of the earlier formers will be more massive. In this 
way the difference in clustering strength is a natural result of hi- 
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erarchical mass assembly in halos and downsizing in galaxies (c.f. 
Foucaud et al. 2010, in press). 

In addition to the relative differences in clustering strength at a 
fixed redshift, we also find a potential difference in how the passive 
and star-forming samples evolve with redshift. The star-forming 
samples follow the behaviour we would expect as a result of the 
downsizing scenario for galaxy evolution. Sub-samples of a given 
K-band luminosity (stellar mass) taken from this population exhibit 
a decline in clustering strength towards 2 = 0. The most lumi- 
nous passive galaxies show tentative evidence for constant cluster- 
ing strength across 0.5 < z < 1.5. A constant value of ro of such 
magnitude indicates that the hosting dark matter halos are no less 
massive at lower redshift. 

Direct comparison with the results of previous work is ex- 
tremely difficult as the UDS is currently unique and this study 
is the first of its type utilising K-band selection. Qualitatively, 
our results agree very well with the findings of Coil et al. (2008), 
iMcCracken et alj fc008h and iMeneux et all fcood) : each of these 
studies finds a longer clustering lengths for red or passive galaxies 
over blue or star-forming galaxies below z = 1. We furthermore 
agree with previous works that any luminosity dependence on clus- 
tering in the sample is minimal in comparison with that of passivity 
over t he range of galaxy pro perties that we have investigated. 

IMcCracken et alj ( 120081) find that at intermediate redshifts 
(0.5 < z < 1), galaxies fit by an early-type template possibly 
have a inverse luminosity dependence, with less luminous galaxies 
having slightly longer correlation lengths. This behaviour is well 
known at low redshift, where passive dwarf galaxies exhibit strong 
clustering and are found to be associated with clusters of galaxies 
dConselice et alj|2003l : IZehavi et alj|2005h . Our measurements for 
red galaxies at z < 1 are not inconsistent with these findings. 

4.1 Discussion of uncertainties 

Throughout this work there are uncertainties that are extremely dif- 
ficult to quantify. Though dust attenuation is taken into account 
during template fitting, there is a well-established degeneracy be- 
tween the stellar age and dust content. Even with the wide range of 
photometry available, we cannot fully overcome this degeneracy. 
In addition, the typical redshift uncertainties are ~ 0.05(1 + z) 
and typical uncertainties in U-B colour are 0.15 mag nitudes. We 
use the ERO definitions o flPozzetti & Mannuccil 1(2000) to examine 
whether our passive selection method is robust. Each of the passive 
objects with best-fit redshifts within the range 1 < z ^ 2 meets 
the i - K > 4 (Vega) ERO criterion, with 1423 of the 1457 objects 
lying in the passive ERO (pERO) region, defined by i,J,K colours. 
In addition, when plotted in the (B — z)ab — (z — K)ab plane, 
our conservative passive sample populates the areas in which we 
wo uld expect to find passive galaxies (see the Appendix, Figure [7] 
and lLane et~ai]|2007h . We are therefore confident that the SED-fit 
star formation properties we use are relatively robust, in compari- 
son with previous broad band techniques. 

Even with robust colours and selection criteria, the photomet- 
ric redshift technique can introduce uncertainties into our measure- 
ments by causing a galaxy to be assigned to the wrong sub-sample 
in redshift. Passive galaxies have a strong break implied by our se- 
lection process and hence the photometric redshift determination 
is more reliable for these galaxies than for star-forming galaxies. 
There is an ongoing ESO Large Programme, the UDSz survey, to 
obtain ~ 4000 redshifts of galaxies at z > 1 in the UDS field. 
Early indications suggest that the photometric redshifts for our pas- 
sive galaxies are indeed highly reliable, with error estimates of 



S z ~ 0.02(1 + z) (though at the time of writing the majority of 
these are at z < 1.75). 

Photometric redshifts are particularly difficult to obtain at 
2 > 1.5. It is precisely beyond this point that the red and blue sam- 
ples converge, due in part to a decrease in the red galaxy clustering. 
On the other hand, the blue and star-forming samples' clustering 
strengths increase over this range. Sub-sample misassignment, par- 
ticularly in redshift, will tend to dilute the clustering signal. For 
the observed convergence to be driven by errors in the photomet- 
ric redshifts would therefore require the clustering signal from our 
red sample to be washed out, while that from the blue sample is 
not. This in turn would require the redshifts of our blue and star- 
forming galaxies to be much more accurate than our red galaxies. 
We suggest that such a case is highly unlikely, but await a substan- 
tial sample of redshifts to confirm this. 

Of greater concern is the possibility that due to photometric 
redshift errors, star-forming galaxies are more likely to be assigned 
to the wrong redshift bins. This would dilute the clustering of the 
blue galaxies and artificially enhance the difference between pas- 
sive and star-forming samples. If dilution due to redshift errors is 
significant we might expect that our clustering measurements for 
star-forming galaxies are smaller than other studies of star-forming 
galaxies at these redshifts. In H08 the clustering strength of sBzK- 
selected galaxies for a sample of median redshift ~ 1.5 cut at 
K = 23 was found to be ~ 6.75 /i~ 1 Mpc. This value is simi- 
lar to the star-forming sub-samples of similar redshifts presented 
here. Although selections in other bands are not directly compara- 
ble we have found that clustering strength is only weakly dependent 
on K-band luminosity. Any moderatel y bright sta r-forming sample 
can therefore be used as a comparison. ICoil et alj (2008) found that 
star-forming galaxies at z — 1 have modest clustering strengths 
of ~ 4ft. - Mpc, with similar values found by Adelberger et al 



2004) (BM/BX-selected galaxies, 2 = 1) and McCrack en et al. 
2008) (2 = 0.6). The z 



Heinis et al 



1, UV-selected galaxies in 
(2007) have a slightly larger correlation length (4.1ft. - x Mpc < 
ro < 5.5ft -1 Mpc), similar to those we report. At similar redshifts 
we find slightly longer (3ft -1 Mpc < ro < 5ft -1 Mpc), but con- 
sistent, correlation lengths compared with most other studies. 

A second method we can use to test this poten - 
tial issue was also used in iMagliocchetti & Maddoxl ( 119990 . 
iMagliocchetti & Maddoxl 119991) noted that there is a systematic er- 
ror introduced into the ro values when it is assumed that the redshift 
distribution is accurately known. They take account of this error by 
assuming the errors in the redshifts of their sub-samples are Gaus- 
sian random and apply a co-efficient, [(12a 2 /Az 2 ) + l] 1 ' 27 , to 
broaden the redshift distribution. In this work modelling the indi- 
vidual n(z) bins as top-hat functions would be inaccurate and so 
we have not attempted to adjust our ro measurements in this way. 
However, this method can still provide an estimate of whether our 
result is driven by redshift errors. Each galaxy has an associated er- 
ror in redshift, derived from the % 2 distribution, which we assume 
here can be represented by a simple Gaussian. In this way each 
galaxy in a sub-sample has a probability of being at a given redshift 
and these individual probabilities are then binned and summed for 
all galaxies within a sub-sample. The resulting distribution is then 
used in place of the original n(z) during deprojection. We find that 
the ro values increase, as expected from a broader redshift distribu- 
tion. However, the difference in clustering between the star-forming 
and passive sub-samples remains, though the difference in implied 
halo masses is slightly reduced. 

This method is very much a simplification and the errors as- 
sociated with the photometric redshifts are likely to be much more 
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complex. In particular the errors on the star-forming galaxies' red- 
shifts may deviate significantly from a Gaussian. Photometric red- 
shift codes are known to favour certain redshifts, introducing ar- 
tificial clumping in redshift space. The spectroscopic redshifts ob- 
tained from the UDSz survey will be vital in understanding possible 
biases arising from using photometric redshifts. We conclude from 
our current tests, however, that errors in the redshift determination 
are unlikely to be responsible for the results that we have found. 

A further, and potentially important source of error, is that of 
cosmic sampling variance. In contrast to H08, we have chosen not 
to take the effects of sampling variance into account. The method 
employed in H08 consisted of splitting the field into 4 and comput- 
ing a measurement for each sub-field. This method provides a esti- 
mate of the sample variance for a field only one quarter of the size 
of the original. The ideal method for estimating the sample vari- 
ance is to repeat these measurements for comparable, independent 
fields. To facilitate such work in the future we re-define our sam- 
ples, using simple photometric selection in the Appendix. However, 
we also note that the clustering of sub-samples at different redshifts 
are largely uncorrelated, yet we find consistent differences between 
passive and star-forming galaxies for subsamples with z < 1.5. We 
suggest that it is unlikely that we would find such consistent dispar- 
ity in the clustering results due to cosmic sampling variance alone. 



5 CONCLUSIONS 

We have used the UDS, the deepest contiguous near infrared survey 
currently available over an area of ~ 0.7 deg 2 to investigate the 
clustering strengths of galaxies over the range < z < 3. Best 
fitting templates, including a dust contribution, have been measured 
for each galaxy. From these templates the rest frame U-B colour, 
absolute K-band magnitude, stellar age and e-folding time of star- 
formation, r, are derived. These quantities are then used to define 
'red' and 'blue' galaxies, a conservative passive sample and active 
star-forming sample. 

Each of these samples are then further sub-divided by red- 
shift and K-band luminosity. Angular correlation function measure- 
ments are made for each of these sub-samples and power laws are 
fit over the scales 9u m < 6 < 0.2 degrees, where 8u m corre- 
sponds to 0.3ft. - 1 Mpc. This limited fitting range is used so as to 
avoid the bias caused by multiple halo occupation (on small scales) 
and to avoid the scales where the correlation function and integral 
constraint correction become similar. Using photometric redshift 
distributions, the angular correlation functions are deprojected to 
find real space correlation lengths. 

We find that the conservatively-defined passive sample is more 
strongly clustered than the star-forming sample for all luminosities 
at z < 1.5. Above this redshift the clustering strengths appear to 
converge, though we would require greater depth and statistics to 
make a firm conclusion. Furthermore, by relaxing our strict pas- 
sive galaxy selection we show that passive galaxies are the most 
strongly clustered of the 'red' galaxies. 

Clustering strength is intimately related to the typical mass 
of dark matter halo hosting a sample, and K-band luminosity is a 
reasonable proxy for the stellar mass of a galaxy. We have therefore 
shown that passive galaxies of a given stellar mass are hosted by 
more massive halos at z < 1.5. At higher redshift the convergence 
of clustering strengths suggests that the red sequence is in the early 
stages of being populated. 

We have also studied how the clustering of passive and star- 
forming sub-samples compares with those of dark matter halos of 



varying mass. The clustering strengths of star- forming galaxies de- 
cline towards z — for a given K-band luminosity (stellar mass), 
and hence are typically found in less massive dark matter halos. 
This finding is consistent with the wealth of evidence supporting 
downsizing. 

The UDS imaging project is ongoing, and expected to gain at 
least 1 magnitude in J, H, K depth by 2012, and potentially substan- 
tially deeper. This will enable us to extend this work to higher red- 
shift and lower luminosities. In addition, our ongoing spectroscopic 
programme (UDSz) will provide several thousand spectra over the 
redshift range probed by this analysis. This will enable a variety 
of more detailed studies, in addition to substantially improving the 
reliability of photometric redshifts. The addition of similar deep 
fields from upcoming surveys (e.g. Ultra VISTA) will also allow ro- 
bust estimates of the sample variance. We therefore anticipate ma- 
jor progress in our understanding of the galaxy populations at the 
crucial epoch when the red sequence first becomes established. 
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APPENDIX 



In this appendix we re-visit the 2-colour selection of iDaddi et al.l 
( 2004). We aim to derive a series of simple boundaries which can be 
used to approximate the passive samples used earlier in this work. 
Our purpose here is to facilitate the reproduction of our results in 
complementary fields rather than create a robust definition for a 
passive sample. We choose the ( B-z')ab-(z ' — K)ab plane be- 
cause the pBzK selection of Daddi et al. (2004) is well understood, 
and also because of the track identified in lLane et alj (2007). This 
track has colours consistent with the SED of an E/S0 type galaxy, 
evolved through redshift, and is well separated from the rest of the 
distribution. It suggests that it may be possible to separate passive 
and star-forming galaxies over all redshifts with relative ease. 

The use of such simple photometric sample definitions may 
be preferable to full star-formation history fitting when comparing 
different data sets. In minimising the complexity of selection, the 
intrinsic differences in the samples should be more apparent and 
hence a better estimate for cosmic sample variance can be obtained. 
The colour boundaries outlined below are defined using the Subaru 
and UKIRT filters. Conversion to other filter sets can be made with 
reference to lHewett et al.l (2006) and Furusawa et al.l (2008). 



Samples in the (B — z) — (z' — K) plane 

The plane is split into cells, 0.2 x 0.2 in colour, and within each 
cell the fraction of galaxies that are in the passive sample is com- 
puted. Boundaries are chosen to provide a simple selection tech- 
nique which includes those cells with > 50% passive galaxies. 
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Figure 7. (B — z')ab — {%' — K)ab colour-colour plot for actively star- 
forming sample (blue points), our 'conservative' passive galaxies (red) and 
those which fall between these criteri a (black). A l so sh own are the com- 
monly used BzK selection criteria of Daddi et al. 12004) and the borders 
defined in the appendix. Only half of the star-forming points are plotted for 
clarity. 



Table 1, Values for the boundary criteria defined in the text. 



Co-eff 


0.l%SFRi 


0.5%SFRi 


l%SFRi 


!0%SFRi 


a 


0.8 


1.0 


1.5 




P 


-1.0 


-1.3 


-2.5 




7 


3.6 






2.2 


S 


-1.4 


-1.2 


-0.8 


-1.0 


e 


7.3 


6.1 


4.4 


4.3 



These borders are shown in Figure[7]together with the conservative 
passive sample and star-forming samples used in the main body of 
this paper. The boundaries for each sample are: 

{z -K) AB < a(B - z')ab + /3 
or (B - z')ab > 7 
or 0' - K) AB > S(B - z') AB + e; (6) 

where the values for a, ft, j, 8 and e are given in tableQ] 

These 2-colour-defined passive samples include all galaxies 
in the relevant region regardless of their best-fit age, U-B colour 
or SFR. The number of objects in each sample are 5438, 7287, 
8986 and 12242 respectively, and the number of objects in the star- 
forming sample is 20645. 

The clustering properties of the these galaxy samples were 
treated in the same way as those in the main body of this 
paper. In each case the properties of galaxies selected by the 
photometrically-defined borders, shown in Figure[7] reflect those of 
the intended population: the majority of sub-samples having very 
similar clustering strengths. Each of the conclusions drawn previ- 
ously are equally applicable to these photometrically-selected sam- 
ples. 

The border-defined passive samples include slightly larger 
numbers of galaxies than the more physically motivated selections 
presented earlier in this paper. The additional galaxy numbers en- 
able measurements beyond z — 2 for some of these samples. 
Though highly uncertain, these measurements indicate equivalence 
of clustering strength at z > 2, as was suggested in section|4] 
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Figure 8. Clustering strengths of the 2-colour-defined galaxies. The symbols have the same meanings as in Figure|4] The values for 
rrj are largely consistent with those in Figure|4] Use of the boundaries defined for these samples should therefore yield measurements 
suitable for comparison with our main results. 
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